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ABSTRACT 


This paper describes aircraft observations made with the NOAA P-3 research air- 
craft on 24 March 1989. The measurement region was over the marginal ice zone, 
southeast of Spitzbergen, between 74 to 76° N and 19 to 27° E. Above the atmospheric 
boundary layer, the geostrophic wind was 12 ms~' at approximately 15° off-ice. A well 
mixed laver extended along the entire 200 km flight path which was perpendicular to and 
centered over the ice edge. There was a stratocumulus layer over the ocean which de- 
creased in thickness towards the ice. The associated inversion, which coincided with the 
cloud layer top, continued to decrease in height over the ice. This caused a strong 
thermal wind effect within the atmospheric boundary layer. The near-surface 
geostrophic wind decreased from approximately 12 ms-! over the ocean to 5 ms“! over 
the ice due to the horizontal temperature gradient and sloping inversion. 
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I. INTRODUCTION 





A. ARCTIC ICE EDGE METEOROLOGY 

Desire to make use of the Arctic polar region resources has greatly increased man’s 
activity there. Understanding the environment is vital in order to tap the immense dis- 
coveries of oil, natural gas and-rich fisheries. The relationship between global weather 
and the polar ice cap has also been an item of long standing interest. In the Arctic, 
between 1953 and 1977, the average maximum extent of the sea ice was 14 million km, 
which is 5%0 of the northern hemisphere. The minimum extent was 7 million km? (Walsh 
and Johnson, 1979). How this annual variation in ice cover is related to large-scale and 
mesoscale circulations is not fully understood. 

The, marginal ice zone (MIZ) is the transition region between the open ocean where 
there ly no sea ice and the polar ice pack with 80% or greater ice cover. Large horizontal 
gradients in ice, ocean and atmospheric properties characterize the region. Also, strong 
interactions between the ice, ocean and atmosphere occur. The spatial variation of the 
atmospheric boundary laver (ABL) thermal structure and surface roughness differences 
between the ice and occan cause secondary atmospheric and oceanic circulations near 
the ice edge. 

Warm ocean currents limit the extent of the ice cover. The wind controls the con-— 
centration of ice in the MIZ (Davidson and Guest, 1987). The ice characteristics from 
the polar ice cap across the MIZ vary. Overland et al. (1983) describe the MIZ by three 
ice regions: 1) a smooth, interior ice pack region of large ice floes with 1 km horizontal 
and 2 to 3 m vertical dimension; 2) a transition zone of small rectangular floes with 20 
m horizontal and 3 to 6 m vertical dimension; and 3) a rough zone of small, broken floes 
with less ice concentration than zone two. These varying conditions of surface 
roughness affect the wind and interactions between the atmosphere, ice and ocean. 

The ice is relatively fresher than the ocean. The melting and freezing of the ice can 
affect the temperature, salinity, and vertical density structure of the upper ocean. The 
change in surface stress across the MIZ is believed to cause ocean upwelling near the ice 
edge (McPhee, 198§3). 

The atmosphere’s temperature markedly varies across relatively short distances, 
from warmer temperatures over the ocean to very cold over the ice. Atmospheric st°- 


bility can also experience abrupt variation. Over the ocean, the atmosphere is typically 











‘unstable, while over the ice it is stable to slightly unstable. Strong inversions are com- 
mon. Shapiro and Fedor (1989) observed very intense cyclones, referred to as polar 
lows, which formed as an arctic cold-air outbreak and cold front moved from the MIZ 
out over the warmer ocean. 

Observations in the MIZ are limited by the region’s remoteness and the resulting 
logistic difficulties that this imposes. Various multidisciplinary experiments have made 
observations of interactions -between the ice, ocean and atmosphere. The Arctic Ice 
Dynamics Joint Experiment (AJDJEX) was a five year project in the 1970s to determine 
ice pack motion and its causes. Using AIDJEX data, Thorndike and Colony (1982) 
studied the causes of average ice motion. They found the long-term average wind and 
the upper-level circulation of the ocean equally contribute to the long-term (several 
months) average ice motion for distances greater than 400 km from a land boundary. 
The ice motion fluctuations (after removing the mean motion) are more closely related 
to the geostrophic wind than to the ocean current. The geostrophic wind accounted for 
10% of the variance of the ice velocity. Within 400 km of land, the effect of the coastline 
on internal ice stresses also contributes to the ice motion. From the Norwegian Remote 
Sensing Experiment in 1979, Johannessen et al. (1983) found a correlation between wind 
direction and ice edge convergence/divergence. With on-ice flow the MIZ is more 
compact than with off-ice flow. 

The most extensive documented ABL observations in the MIZ have come from the 
Marginal Ice Zone Expe,siment (MIZEX) in the Bering and Greenland Seas during 1983 
and 1984. Johannessen (1987) describes the summer MIZEX East program in Fram 
Strait and the Greenland Sea. Fairall and Markson (1987) show the variation of surface 
fluxes for five aircraft flights during July 1983- over the MIZ in the Greenland Sea. 
Conditions for the five flights include winds from all four directions relative to the ice 
edge at Wind speeds 7 ms~! and below. They noted that the wind direction and its effect 
On atmospheric stability, can be an important factor in describing the surface stress. 
On-ice flow or flow parallel to the ice edge from the south at less than 5 ms“ , caused 
more stable conditions in the MIZ. Under these conditions, the effects of hydrostatic 
stability on surface stress were the same order of magnitude as the effect due to varying 
surface roughness. Kellner et al. (1987) document turbulence statistics and spectra un- 
der conditions of off-ice flow in the Fram Strait. This study includes aircraft and 
shipboard observations made during the summer of 198-1. The mixed layer depth in- 
creased towards the ocean as convective mixing from the warmer occan became impor- 
tant. This coincided with a shift in the peak frequency of the vertical wind component's 


variance spectra from | Hz (100 m) over the ice to .4 Hz (250 m) over the ocean. Sec- 
“ondary low frequency peaks were evident over the ocean. These indicated secondary 
circulations. 

Early numerical models of the MIZ simulated off-ice flow. Overland et al. (1983) 
used a one-layer slab model to show the evolution of the ABL. The wind velocity in the 
ABL and the height of the base of the inversion proved sensitive to the surface 
roughness decrease and heat flux increase from ice to ocean. By 50 km oceanward of 
the ice edge, the horizontal temperature gradient had caused an 8% increase in wind 
speed. Roughness changes were responsible for a 9% wind speed increase across the 
MIZ (approximately 35 km). The base of the inversion slowly increased towards the 
ocean. This study agreed with the observation that an off-ice wind ficld contributes to 
ice edge divergence. Reynolds (1984) enhanced the slab model by including cloud and 
radiation effects. The addition of cloud and radiation effects for off-ice flow, slightly 
increased the slope of the inversion base and slightly decreased the wind speed. Such a 
slab model adequately simulates the shallow boundarv laver of off-ice flow, but for flow 
from other directions, a more sophisticated model is needed. 

Benuctt and Jfunkins (1986) used a multilevel, two-dimensional, higher-order closure 
model to simulate on-ice flow. The model included stratus cloud effects and sea ice 
thermodynamics. They used cuserved data to set the model’s initial conditions of wind 
and potential ternperature profiles at the ice edge. The processes of turbulence, lateral 
and vertical advection and radiation were simulaied. The model showed good agreement 
with the observed vertical soundings at various distances from the ice edge. The base 
of the inversion increased from 500 m at the ice edge to 1000 m at 150 km iceward of the 
icc edge. The study confirmed the importance of including cloud top longwave cooling 
when modeling a cloud topped ABL under conditions of on-ice flow. The model results 
suggested that the boundary layer cooling away from the ice edge was mainly duc to 
cloud top longwave cooling, not downward heat flux. 

Kantha and Mellor (1989) used a two-dimensional, higher-order closure model to 
simulate off-ice, on-ice, and along-ice flow. For off-ice flow, the model showed good 
agreement with the Overland et al. (1983) results. With on-ice flow where there was 
strong heat flux from the ABL to the ice, an inversion tens of meters in depth and a 
strong jet resulted. Yor along-ice flow with a slight component either on or off-ice, a 
strong thermal ice edge front with intense turbulence resulted. The number of detailed 
descriptions of the ABL over the MIZ from observations are few. More observations 
are needed te validate the model results. 


B. CEAREX 

The Eastern Arctic Aircraft Expedition, which took place in March 1989, was a 
component of the multinational, interdisciplinary Coordinated Eastern Arctic Exper- 
iment (CEAREA). This expedition provided the opportunity to investigate the structure 
of the ABL in the vicinity of the MIZ in the Barents Sea using the NOAA P-3 research 
aircraft. 

Under conditions of flow nearly parallel to the ice edge, one flight was made on 24 
March 1989 to make measurements of the ABL, the background atmospheczic structure, 
and turbulent fluxes using a gust-probe. This paper discusses the observations of that 
flight. This is the first set of observations to provide a comprehensive description of 
ABL structure over the MIZ for flow nearly parallel to the ice edge. 





Il. THE DATA 
A. FLIGHT DESCRIPTION AND CONDITIONS 


The U.S. Nasy Fleet Numerical Oceanography Center’s surface pressure analysis 
of 1200 Coordinated Universal Time (UTC), 24 March 1989 (Figure 1) depicted two low 
pressure centers, one west of Norway at 60° N and 0° E and one at approximately 73° 
N and 50° E. There was a short inverted ridge where the isobars and wind trajectories 
curve south between these two low pressure centers. In the vicinity of the short inverted 
ridge, which included the measurement area centered at approximately 75° N and 23° 
E, the wind direction varied. The flov’ in the measurement area was northeasterly. To 
the cast of the measurement area the flow was more northerly, while to the west it was 
easterly. 

During the aircraft's approach to the measurement area, stratocumulus clouds were 
observed over the open ocean wih some precipitation beneath the clouds. This cloud 
laver thinned towards the ice edge and there were no clouds over the Ice. 

The criterion for the flight was that the wind flow would be nearly parallel to the ice 
edge, with a small off-ice component so boundan layer operations would not be hmited 
by fog. With the criterion met, measurements were made in the MIZ southeast of 
Spitzbergen on 24 March 1989 berween 74 to 76° N and 19 to 27° E. The MIZ itself 
was fairl. compact. The ice pack (with greater than 80% ocean cover) followed a 19 to 
20 km zone of pancake ice. The ice edge was oricnied northeast-southwest as indicated 
by the dashed line in Figure 2. 

The aircraft speed during measurements was approximately 100 ms"! . The flight 


pattern (Figure 3) included: 


¢ Four 200 km Iegs perpendicular to the ice edge at 900, 800, 700 2b and 45 m which 
were almost exactly centered over the MIZ. 


¢ Two legs approximately 30 km in length over the ice at 33 and 130 m, parallel to 
the ice edge. 


¢ Two legs approximately 50 km in Ienath over the ice at 73 and 150 m, perpendicular 
to the ice edge. 


tH 
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1200 UTC 24 March 


Figure 1. Surface pressure analysis 
Source: Lackmann et al., (1989) 


e Two Icgs approximately 50 km in length over the ocean at 45 and 267 m, parallel 
to the ice edge. 


¢ ‘Two legs approximately 50 km in length over the ocean at 135 and 267 m, per- 
pendicular to the ice edge. 


B. INSTRUMENTATION 

The ABL mean and turbulent variabics were-measured from the NOAA P-3. Ob- 
taining a wind measurement with respect to the earth requires determining both the 
aircraft's velocity with respect to the carth ait the velocity of the air with respect to the 
aircraft. An inertial navigation systein (INS) provided the aircraft velocity and position 
with respect to the earth. The mean velocity of the air with respect to the aircraft was 


obtained using a pitot-static tube. A forward port on the tube gives a total pressure 
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Figure 2, NOAA-10 IR imagery of 1038 UTC 24 March 1989 


measurement and side ports give static pressure. The difference between these two 


pressures is a measure of dynamic pressure. Combining this dynamic pressure meas- 
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Figure 3. Flight pattern schematic for CEAREX flight 


urement with ambient temperature and static pressure measurements gives the mean 
wind speed. The wind direction was obtained from vanes attached to the tube. The 
turbulent wind measurements were obtained using a gust probe whose performance is 
described by Greenhut and Gilmer (1985). 

A. Rosemount votal temperature probe was used to obtain the mean temperature. 
Total temperature refers to the ambicnt temperature plus dynamic heating (the effect 
of air flow on the sensor). A .025 mum platinum wire sensing clement experiences 


changes in specific resistance which are proportional to the change in temperature. This 


same principal is used by the bead thermistor which gave the temperature measurement 
as part of the turbulence package. The sea surface temperature was remotely sensed by 
a Barnes PRT-5 radiometer. It compared the radiance of the sea surface in the 9.5 to 
11.5 wm wavelength band with that of an internal reference cavity. 

A dewpoint hygrometer provided an absolute measure of humidity (mass per unit 
volume). A mirror is maintained at the temperature where the moisture in the air just 
begins to condense on it. The condensation is monitored optically. This sensor was 
slow to respond to a change in the measured variable compared to the other sensors. 
It took tens of seconds to reach equilibrium. Another part of the gust probe system, a 
microwave refractometer, measured the turbulent moisture fluctuations. This sensor 
directly measures fluctuations in the radio refractive index. Combining this measure- 
ment with those of temperature and static pressure gives the water vapor density. 

There was also a lidar system on board the NOAA P-3 which was operated by 
University of “Vashington personnel. The system transmitted two laser light beams ev- 
ery second, of wavelengths .532 and 1.046 um. The transmitted light is strongly scat- 
tered by cloud droplets. A tclescope collects the returned light and focuses it on a 
detector. From there it is amplified, digitized every 20 ns and recorded. This resulted 
in a resolution of 7.5 m. Cloud base and cloud top height measurements were obtained 


using the lidar. Table | lists the variables, their sampling rates and the sensors used. 


Table 1. INSTRUMENTATION CARRIED BY THE NOAA P-3 IN CEAREX. 


Mean Variables Sampling Rate (s~!) 


: 


bined with INS 
Dewpoint 












Rosemount Total Tem- 
perature Probe 
1O2CH2AF 


| 1 General Eastern 120DEP 
Dewpoint Hygrometer 


Altitude J Radar altimeter (Stewart- 
Warner APN159) 
Pressure i Rosemount 542K2 


Barnes PRT-5 
Radiometer 
Turbulent Variables Sampling Rate (s~') Instrument Source 
u, V, W 40 Gust probe combined 
with INS 
i 
q 40 
















Sea surface temperature 






C. DATA PROCESSING 
Most of the mean variable sensors were checked by comparing their measuiements 
against routine observations taken at Bodo, Norway, where the CEAREX flights origi- 


nated. This comparison revealed: 
¢ The dewpoint hvgrometer observations were 5.6° too cold. 
e Temperature data were within 1° K and pressure data within 1 mb. 


The radiometer uscd to obtain the surface temperature was calibrated at NOAA‘s Office 
of Aircraft Operations (OAO) in Miami, prior to and after the CEAREX flights. These 


checks showed a drift of 2° K at -35° C and I[° K near freezing. 














Inertial navigation systems (INS) generally experience slowly varying drifts with 
time. It was planned to use the Loran-C navigational aid signal to compensate for this 
following the flight. However, tne Loran-C signal was too weak to use. So, the possi- 
bility exists that INS drift could have introduced variations as large as 3 ms=! in the 
wind field (Shaw, 1988). 

Greenhut and Gilmer (1985) described the calibration and performance of the gust 
probe system. The thermistor is calibrated in a wind tunnel to correct for dynamic 
heating. This correction is approximately 4° C at an air speed of 100 ms=', The char- 
acteristics of the microwave refractometer are also well known. The gust prebe system 
performance during a set of aircraft maneuvers (called the Lenschow maneuvers) is well 
established. During the CEAREX flight, these Lenschow mancvvers were executed. 
Later evaluation of the gust probe system data during these maneuvers identified any 
abnormalities in the performance of the system. 

The gust probe system was mounted on a boom attached to the nose of the aircraft. 
The boom length was chosen so its natural frequency is about 14 Hz. Freqnencies 
greater than 11.5 Hz are filtered out of all the gust probe system signals on-board the 
aircraft, prior to recording the data. This avoids aliasing the higher frequenc: natural 
resonances of the boom and other sensors into the recorded data. 

Two separate data streams of mean and tuibulent (gust probe) observations were 
recorded, with UTC being the only common variable. The mean data were initially 
processed at NOAA’s Office of Aircraft Operations (OAO) in Miami. All the above 
noted errors were corrected and the data were screened for unrealistic values. 

The turbulence data were initially processed at NOAA’s Environmental Research 
Laboratories (ERL) in Boulder, Colorado. The turbulent fluxes were determined using 
the eddy correlation method at the Naval Postgraduate School. It was desired to show 


the variation of the boundary laver fluxes compared to the surface variability and also, 


I] 


to use Fast Fourier Trensforms (FFT) which use less computer time than other methods. 
To accomplish this, the data (time series) were analyzed in 102.4 s segments corre- 
sponding to about 10 km of sampling distance: This 102.4 s segment length met the FFT 
requirement of 2’ data points per segment. The sampling rate was 40 s-', yielding 4096 
( 2" ) points per Fourier transform. 

As noted earlier, the recorded gust probe system data included frequencies. up to 11.5 
Hz. The first step was to pass the time series through a high-pass symmetric compound 
running mean subtraction filter (Holland, 1968). This eliminated the unwanted low fre- 
quency (mean) information from the time series, leaving the high frequency (perturba- 
tion or turbulent) part. The filter parameters were set such that more than 90% of signal 
variance was removed for wavelengths longer than 10 km (Shaw et al., 1990). This is 
equivalent to the filter eliminating 90% of the signal information for frequencies less 
than .0] Hz. This particular filter was chosen for its ef“ictency anu ease of implementa- 
tion. Another benefit of using this particular filter vas that its transfer function is 
well-behaved and does not seriously distort the other frequencies i:: the time series. The 
filter passed a small amount of signal variance for frequencies of .01 to .03 Hz, and all 
signal variance for frequencies greater than .03 Hz. It was desired to retain signal in- 
formation for frequencics of .01 Hz and greater. The information that was eliminated 
for frequencies .01 to .03 Hz was reinserted during a later step. 

Next, Fast Fourier Transforms (FFT) were performed tn tach 192.4 s segment and 
the spectra and cospectra were calculated. The one-sided spectra and cospectra, 
G..(/) and G,,(/) where a and b are are any two variabies and fis frequency, describe the 
Variance and covariance distributed cver frequency. The variance and covariance are 


defined: 


i= | Gada 


a’b! = | Gafaf 
0 





The covariances represent turbulent fluxes (Sovbjary 1989). For example; the vertical 


temperature flux (w’T’) is a covariance between tic .” ulent or perturbation quantities 





«(denoted-by the prime) of vertical motion (w’) and tengerature (T’). 

The third step was to divide the spectra and cospectra by the transfer function of the 
compound running mean subtraction filter, which was apflied'in the-time domain. The 
transfer function describes the relation between the Signa! input to the filter-and the fil- 
ter’s output. It is obtained by performing a Feurier Transform onthe function de- 

; scribing the filter. The covariance transfer function is a ratio between the covariance 
after filtering and the covariance before filtering. Tne cospectra computed in step two 
, is a ratio between co-ariance after filtering and frequency. Thus dividing the cospectra 
from step two by the transfe: functicn gives the true cospectra (cospectra before filter- 
ing). This relation is also valid for the spectra if variance is substituted for covariance. 
This step reinserts the .01 to .03 Hz information that was removed by the filtering instep 
one. 
The last step was to integrate the true spectra and cospectra with respect to ffre- 
quency. This gave the desired average variances and covariances (turbulent fluxes) for 


each 102.4 s segment. 
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UI. MEAN STRUCTURE 


A. VERTICAL STRUCTURE 

Vertical profiles were obtained from-iU's averages during altitude changes between 
the straight and leve: measurement legs. These altitude changes took place at apprvxi- 
mately 100 km from the ice edge over.the ice a1.d ocean. The composite virtual potential 
ternperature prcfile (Figure 4) showed that the ABL was well-mixed at 2 yproximutely 
252° K to about 150 m over the ice and at approximately 263° K to abr ut 750 m-cver 
the ovean. Above the ABL, the profiles over both the ice and ocean display w:early 
identi-al thermal structure (+2° K). he composite specific humidity (q) profile 
(Figure 4) had this same general structure except that over the ocean, the humidity 
profile was not so well-mi. ed in the boundary layer. 

Four altitude changes were made over the ice. The resulting profiles individually 


revealed: 


e Three of the four virtual potential temperature profiles were well-mixed to about 
100 m (Figure 5\. There was a slightly stable layer between 100 and 200 m, ana « 
strong capping inversion above 200 m. Ths fourth profile which was the farthest 
from the ice edge differed in that the well-mixed layer was approximately 70 m deep 
and inversion layer was .5° K warmer than the other thrce profiles. This suggests 
the inversion was sloping downward sharply toward the ice. 


e The de «point hvgrometer had a tendency to overshoot a step change. Two of the 
four available specific humidity (q) profiles (Figure 6) were obtained from altitude 
changes that started within the ABi, Tiese showed a constant q value, charac- 
teristic of a well-mixed layer, to a height of 150 m. The other two profiles, one of 
which was the profile farthest from the ice edge, were obtained from altitude 
changes that started above the ABL. These showed-q variation within the ABL, 
making identification of the well-mixed layer top less clear. However, for the pro- 
file farthest from the ice edge, q began to decrease with height at approximately 70 
m. This agrees with the virtual potential temperature finding that the ABL depth 
for the profile farthest from the ice edge is about one-half that found in the other 
profiles, 
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Fk sure 4. Composite profiles of a) virtual potential temperature and b) specific hu- 
midity 


e The downstream profiles of virtual potential temperature and specific humidity re- 
flected moistening of the boundary layer and warming due to surface heat flux. 


¢ The wind speed profiles were seriously contaminated by turns and rapid altitude 
changes. Using only 10 s data segments where the aircraft heading changed less 
then .5°s"! and the altitude changed less 3 ms- , the profiles (Figure 7) showed 
thut the ABL wind adjustment to gcostrophic value occured almost entirely in the 
inversion. Between 150 and 250 m, the wind increased by § ms". 


B. CROSS-SECTIONS 
Shaw et al., (1990) described the structure of the ABL during the flight. The data 


from the 200 km flight legs provided cross-sections perpendicular to the ice edge: The 
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Figure 5. Prof‘les of virtual potential temperature over ice 


radiometric surface temperature cross-section (Figure 8) had a step-like structure. Over 
the ocean, the surface temperature was slightly above 0° C at distances greater than 40 
km from the ice edge and just below 0° C within 40 km of the ice edge. The ice edge 
was clearly identified by a marked decline in surface temperature, with leads indicated 
by warm temperature spikes. Twenty kilometers iceward of the ice edge, the radiometric 
surface temperature had dropped to -15° C and by 60 km to near -25° C. 

Virtual potential temperature traces (Figure $) of the 900, 800, and 700 mb Icgs re- 


flected near horizontal temperature contours which indicate there was virtually no 
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Figure 6. Profiles of specific humidity over ice 


change of temperature perpendicular to the ice edge at any altitude above the ABL. 
Therefore, interactions with the stisace are responsible for all cross-ice-edge thermal 
variability within the ABL. Ocean-ice temperature differences cause baroclinity within 
the ABL. Consistent with the small off-ice component to the flow, the baroclinity was 
strongest in the region 40 km oceanward of the ice edge. 

Using the lidar pointing downward on the 700 mb leg and upward on the 45 m Ieg, 
it was possible to determine the top and bottom of the stratocumulus layer. The 


stratocumulus layer (shaded region on Figure 8) was 300 to 400 m thick at 100 km 
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Figure 7. Profiles of wind speed over ice 


oceanward of the ice edge. The cloud base was at approximately 225 m. In the vicinity 
of showers, the cloud buse was lower. The stratocumulus layer gradually thinned to no 
clouds over the ice. The Jidar-derived cloud layer base was comparec to the calculated 
lifting condensation level (LCL). The dashed Jine on Figure 8 identifies the LCL. The 
results were generally consistent, though the LCL was approximately 25 m lower. The 
cloud top was used as an estimate of the boundary layer height (h). This method could 
not be used to determine h over the ice, since there were no clouds. The potential tem- 


perature profiles were well-mixed in the ABL and virtually idenucal above the ABL over 
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the ice and ocean. The boundary layer height (h) was estimated by assuming that the 
45 m leg potential temperature data represented the ABL well-mixed value. Then the 
height where this same potential temperature value was found on the free-atmosphere 
profile was considered the ABL depth. This method assumes a barotropic atmosphere 
above the inversion in the cross-ice edge direction and therefore does not allow for a 
thermal wind component parallel to the ice edge. This assumption is not quite correct, 
since the influence of subsidence from the large scale inverted ridge would act to increase 
the temperature gradient in the sloping inversion and make a thermal wind component 
parallel to the ice edge likely. However, the parameterization appeared to give a rea- 
sonable estimate of h over the ice. The lidar cloud top height observations were in close 
agreement with the parameterized h values near the ice edge. 

The specific humidity traces (Figure 9) showed that within the ABL, the air over the 
ocean was twice as moist (1.2 g kg-') as the air over the ice (.6 g kg-'). The atmosphere 
above the ABL was quite dry with specific humidity values of .4 g kg-! just above the 
ABL and .2 gkg-! at approximately 1400 m. In agreement with the thermal structure, 
the strongest horizontal moisture gradient was in the region 40 km oceanward of the ice 
edge. 

Windspeed (Figure 10) varied by only | ms-' along the 200 km flight legs above the 
ABL. Within the boundary layer over the ocean, the windspeed was within 1 ms-! of 
the windspeed above the ABL. Over the ice, however, the windspeed near the surface 
varied from 11 ms-! at the ice edge to about 6 ms~ at 100 kilometers from the ice edge. 
The strongest horizontal windspecd variation occured in the region more than 40 km 
iceward of the ice edge. 

The wind direction cross-section (Figure 11) showed a variation of approximately 
10° at most, on all horizontal flight legs. The near surface wind was from approximately 


25°, Between the surface and 900 mb (about 800 m), the wind veered no more than 10 
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to 12°. Between altitudes of 1 and 3 km, the wind backed 20°. This implies cold 
advection. However, some of the variation may be due to unaccounted for drifts in the 
INS system. Therefore, the wind direction changes, particularily above the ABL where 


the changes are small, should be regarded with some caution. 


C. HORIZONTAL GRADIENTS 

The flight pattern with its parallel and perpendicular ice edge legs allowed horizontal 
gradients to be determined at two levels over both the ice and ocean. Figure 12 shows 
the spatial series for u, v and virtual temperature and their linear best fit from the 50 km 
segment at the end of the 45 m flight leg over the ocean. u and y are in ice edge coor- 
dinates with positive u towards the ocean. A linear best fit was obtained for each vari- 
able. This provided a measure of how much a variable changed (its delta. value) along 
the flight leg. The delta values were used to determine the components of the horizontal 
gradient using a icast squares method for heights of 40 and 150 m over the ice and 45 
and 267 m over the ocean. Results are listed in Table 2. Using temperature T as an 


example, values of AT, Ax, and Ay derived from each leg were used in: 


ATz ee 


Ax + Ce. Ay 
Ox Cy 


This is of the form z = ax + by. The least squares method minimizes )((z — z,? with 
f 


respect to a and b. 


o 
aa ) te — (ax; + by)? = (0) 
; 


and 
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Dt — (ax; + by) = (0 : 
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Expanding.and rearranging gives: 


and 


_ mie Dye Daas 


The gradient calculation for pressure was slightly different. Figure 13 shows the 
pressure series reflected-the variation of the aircraft’s altitude. Using the radar altimeter 


height (z) values and the hydrostatic approximation, the pressure was adjusted to a 


constant flight altitude (z,,,) then used in the gradient calculation. 
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P= pRTy 


—_ P ‘ 
P RTy ' 





where p is pressure, p is density, g is gravity, R is the specific gas constant, and T, is 
virtual temperature which is assumed constant over small altitude changes. Substituting 


this into the hydrostatic approximation and integrating gives: 





ie ean Ete Oe 
n( Prep ) ae (2 — Z;e)) 
ee ee ee 
Pref exp| RTy (2 Zref) | 


§ 
Pre =P exp] Rly. (z— if | 


where p,,, is the pressure adjusted to a constant flight altitude. - 
Figure 14 is tke specific humidity series from the dewpointer and the microwave 

refractometer (10 s averages). The mean has been subtracted from the refractometer 

data. This is why the mean values differ between the two series. The dewpointer’s 

tendency to overshoot a step change affected the delta values. For this reason, the 

microwave refractometer data were used to calculate the horizontal specific humidity 

gradients. Table 2 summarizes the horizontal gradients over the ice and ocean. The 

gradients are with respect to the ice edge coordinate system where positive x is towards 

the ocean and positive y is along the ice edge towards the northeast. The geostrophic 

Wind was calculated using the pressure gradient. The other gradients were used in de- 


termining the physical balances described in the next chapter. . 
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A ERERSOONE 2. GRADIENTS DERIVED USING THE LEAST SQUARES METHOD 


Over ice Over ocean 
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D. GEOSTROPHIC WIND 




































= (m / S km) 





The geostrophic wind (J’,) was calculated using the gradients from the pressure 


sensor data. 
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Figure 15 shows the geostrophic wind over the ice and occan. iv at 40, 45, 150 and 267 
m was calculated using the pressure gradients. Since 900 mb was at all times above the 
ABL, the actual wind at that level is considered geostrophic. Over the ocean, the 
geostrophic wind decreased in magnitude between 45 and 267 m. This is consistent with 
the ABL thermal gradients. The thermal wind (17) defined by the gcostrophic wind 
change with height, blows with low temperature to the left. Over the ice, the geosti ophic 
Wind changes in magnitude and direcuon between 40 and 150 im. ‘Ths is not consistent 
with the ABL thermal gradients since V, defined by the change of I’ with height does 
not blow with low temperature to the Ieft. Jt is the inversion temperature giadicnt noted 
in Figure 4 that accounts for the I’ change. Above the ABL, I’ over the ice and ocean 
has comparable magnitude and direcuon. Thus, over the ice, I" is relatively constant 
with height within the ABL but changes significantly in the inversion, while over the 
ocean, I’ changes significantly in the ABL. This relation has implications concerning 


the momentum balance described in the next chapter. 
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IV. TURBULENCE MEASUREMENTS 


This section describes the heat, moisture and momentum balances in the ABL ap- 
plied at the corners of the flight pattern over the ice and ocean. An ice edge coordinate 
system is used with positive x towards the ocean and positive y along the ice edge to- 


wards the northeast. 


A. HEAT BALANCE 


The time rate of change of potential temperature, neglecting radiation, 1s: 





where @, u and v are the mean potential temperature and velocity, and #” and @" are the 
turbulent vertical velocity and potential temperature. The coordinate system 1s such that 
X is perpendicular to the ice edge and positive toward the ocean. The advection terms 
represent cross-ice (u o and along-ice advection (+ 2 The time rate of change of 
potential temperature was determined from two flight on along the same course but 
at different altitudes, which were separated by approximately I hr. This is appropriate 
since both flight legs were in the ABL, and the potential temperature profile was well 
mixed. The velocity components, u and 4, were determined from 20 km segments at the 
flight pattern corners. The average vertical turbulent heat flux (::"0") was determined 
(see Chapter II, sections B and C) for the 40 km segments closest to the flight pattern 
corners. Figure 16 shows the s:’0" profiles and linear best fit which was forced to zero 
at .Sh where h is the height of the ABL (120 m over ice, 600 m over ocean). The best 


fit was forced to zero at .Sh since this is the expected height where the heat flux changcs 


sign due to entrainment of higher temperature air from above the version base. 
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Figure 16. Vertical turbulent heat flux 


Table 3 shows the heat balance over the ice and ocean. Over the ice, there was no 
detectavle change in potential temperature with time, while over the ocean, there was 
net cooling. Advection mostly balanced the vertical heat flux. Along-ice-edge advection 
was approximately twice as large as cross-ice advection. The magnitude of the terms 
was larger over the ocean than ice. Compared to the other terms, the residual was small. 


This strongly suggests that neglected effects did not play a large role in the balance. 
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Table 3. HEAT BALANCE 


residual 
(K / hr) 





Source: Shaw et al., (1990) 


B. MOISTURE BALANCE 


The time rate of change of moisture is: 





Gq 0g Cg = Ow'q' 
— = - ue I 
Cl Ox ey Oz 


where q is the mean specific humidity, and q’ is the fluctuating part. The time rate of 
change, u, v and the average vertical moisture flux (w’g’) were obtained by the same 
method used with the heat balance. Figure 17 shows the w’g’ profiles and linear best 
fit which was forced to zero at h (150 m over ice, 750 m over ocean). The best fit was 
forced to zero at h since upward moisture flux is expected to decrease with height up to 
the inversion base. 

Table 4 shows the moisture balance over ice and ocean. The time rate of change 
reflects net drying of the ABL. Over the ocean, advection partially balances the vertical 
moisture flux. The magnitude of the residual indicates this description of the moisture 
balanve was not as successful as with the heat balance. The Figure 17 best fit over the 
ocean does not reflect the observed increase of the moisture flux with height over the 


lower ABL. This increase with height indicates drying near the surface and moistening 


of higher levels. This suggests there is a moisture source in the ABL other than the flux 
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Figure 17, Vertical turbulent moisture flux 


from the ocean surface. A similar increase of w’q’ with height over the ocean was noted 
in a study by Nicholls and Readings (1979). Their observations showed w’q’ increased 
with height over the lower ABL and decreased with height in the upper ABL. The mean 
specific humidity decreased with height throughout the mixed layer. They concluded the 
increasing w"g’ profile was acting to change the decreasing mean profile to a mor well 
mixed (constant) profile. The situation during the CEAREX flight differed slightly in 


that the mean specific humidity profile didn’t decrease with height throughout the mixed 
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layer. However, Figure 4 shows it was not so well mixed at the lower levels. No ob- 
servations were available to determine if w’g’ decreased with height in the upper ABL. 
The moisture flux increase with height implies a downward transfer of water vapor 
in the boundary layer. For this to occur, a moisture source above the surface is required. 
Donelan and Miyake (1973) suggested evaporation of water droplets relea.ed from the 
surface by whitecapping as a possible moisture source that could explain the flux in- 


crease with height. 


Table 4. MOISTURE BALANCE 











residual 
(g / kg hr) 









(g/kghr) | (g/kghr) | (g/kghr) | (g/kg hr) 


Ice -.001709 -.00612 -.001235 -.00788 01353 
“14799 | _.29124 


C. MOMENTUM BALANCE 


The momentum balance, neglecting accelerations is: 





Ouw og 
Gz = flv ty) (J) 
— SE = flu- 1) 2) 


The above equations say that the difference between the pressure gradient force (/v,) and 
the coriolis force (f) is balanced by the momentum flux divergence (stress). Therefore, 
the difference between the observed and geostrophic wind is related to the stress profiles. 
If the difference is constant, the stress changes linearly with height. Ifthe difference isn’t 
constant, curvature in the stress profile results. The observed wind averages for the 


flight levels within the ABL are listed in Table 5. 
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Table 5. OBSERVED WIND AVERAGES 


= 
615 
12.22 
[Ocean 4am SSOS~S dt 
















1. Over the Ice 
Table 5 indicates that u varies and v is constant with height over the ice. Recall 
that over the ice, V. did not vary with height within the ABL since the boundary layer 
was too shallow for the thermal wind to have a significant effect. So, when applying the 
momentum balance in the ABL, V, is assumed to be constant with height. Since v, and 
y are considered constant, their difference is constant and equation (1) requires the 
profile of u’w’ to be linear. Figure 18 shows the ww’ profile and linear best fit. The 
coordinate system is oriented such that the uw’ profile is perpendicular to the ice edge. 
The u wind component varies in the boundary layer. Assuming that u varies 


linearly with height allows u to be described by: 
u=U,+ Su Z (3) 


Where u, is u at the surface and z is the height above the surface. Integrating equation 


(2) over height gives: 


vs0"(z) — v’v'(0) = | ne, — u)dz (4) 


38 





300 


ICE 


FLIGHT LEG 


1 = PERPENDICULAR 
O = PARALLEL 


TO IGE EOGE 


HETGHT (M) 
100 150 200 250 


20 


© 


-0.05 0.05 0.15 


U MOMENTUM FLUX 








0.25 -0,05 0.05 Q.15 0.25 





300 









ICE 





-FLIGHT LEG 
O = PERPENDICULAR 
O ~ PARALLEL 

TO IGE EDGE 


200 


200 


150 


100 


20 


ae) 


V MOMENTUM FLUX (MS) 


. ° % : 
Figure 18. Vertical turbulent momentum flux over the ice 


Substituting equation (3).in equation (4): 


v'se"(2) — v'v"(0) = | Le ~ iO=— > Au Ls la 


Integrating and rearranging vields: 


v'w'(z) = v'w"(0) + fil — u(0) — > 


Ait 


Az 
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(5) 


(6) 


2 


oe, Soe 








Equation (6) describes the best fit curve. Using the v’w’ averages from the various alti- 
tudes, a least squares method was used with this curve to obtain v’w'(0). Table 5 values 
were used for ae , and equation (3) was solved for u(0). Figure 18 shows the vw 
profile and best fit defined by equation (6). The coordinate system is oriented such that 
the vw" profile is parallel to the ice edge. The curve closely matches the observed v'w’ 
averages of the flight-legs perpendicular to the ice edge while there’s some scatter for the 
flight legs parallel to the ice edge. The closeness of the best fits to the observed values 
confirms that the neglected inertia_terms did not have a dominant role. 
2. Over the Ocean 

The large heat flux and surface stress produced a strong mixed layer with little 

change of observed winds with height (Table 5). Recall that over the ocean, the mag- 


nitude of V, varied with height in the ABL (Figure 15) due to the temperature gradient. 


Table 6 hsts the geostrophic wind components over the ocean. 


Table 6. GEOSTROPHIC WIND AVERAGES 


Ocean 44m 






Table 6 indicates that u, is relatively constant and 1, varies with height. Since u, and u 


are considered constant, their difference is constant and equation (2) requires the profile 
of v’w' to be linear. Figure 19 shows the v's’ profile and linear best fit. 


Assuming that v, varies linearly with height allows v, to be described by: 


& 
Y= ‘sot 2 (7) 


where y,, is v, at the surface. Integrating equation (1) over height gives: 
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Figure 19. Vertical turbulent momentuin flux over the ocean 


Ww(2) — ww'(0) = ec — yds (8) 


Substituting equation (7) in equation (8): 


a —— z Ay 
u’sv'(z) — u’w'(0) = \ i) — ¥(0) — ae |e (9) 


4] 





ectoaes Stan Av 
u'w'(z) = u'w'(0) + fl — v(0) — + g | (10) | 


—-2 
Az 


Equation (10) describes. the expected functional form of the variation of u’w’ with-height. 
Table 6 values were used for a , and equation (7) was solved for v,(0). Using the 
u'w’ averages from the various altitudes, a least squares method was used with this curve 
to obtain u’w'(0). Then u’w'(0) was used in equation (10) to obtain the best fit. 
Figure 19 shows the u’w’ profile and best fit. The negative ww’ of the 45 m flight leg 


remains unexplained. There were no indications of equipment malfunction. 
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V. SUMMARY 


Data collected with the NOAA P-3 research aircraft was used to examine the mean 
structure and physical balances of the ABL in the MIZ. The wind flow was nearly par- 
allel to the ice edge with the 900 mb geostrophic wind approximately 15° off-ice. The 
atmospheric conditions above the ABL, were virtually identical at 100 km from the ice 
edge over both the ice and ocean. Therefore, the structure of the ABL was due to 
interaction with the surface. The well mixed boundary layer extended at least 100 km 
into the ice. The horizontal gradients of temperature and moisture were maximum at 
approximately 40 km oceanward of the ice edge but no boundary layer front was evident. 

The heat balance showed that the vertical heat flux was mostly balanced by 
advection. Over the ocean, the magnitude of both the heat flux and advection was quite 
large. The residuals were small compared to the other terms. This strongly suggests that 
neglected effects did not play a major role. 

The moisture balance gave an unusual result. Over the ocean, the profile of vertical 
moisture flux increased with height over the lower ABL. This profile could be viewed 
as acting to make the mean specific humidity profile more constant or it could be de- 
scribing the presence of a counter-gradient flux. 

The relationship between the observed and geostrophic winds in the momentum 
balance explained the curvature in the stress profiles. The boundary layer was suffi- 
ciently deep over,the ocean, to allow the horizontal temperature gradient to affect the 
geostrophic wind. Over the ice the ABL was shallow, thus minimizing this effect. 

These observations provide evidence of the interactions occuring in the MIZ under 
conditions of flow nearly parallel to the ice edge. The sloping inversion and horizontal 


temperature gradient had a significant effect on the geostrophic wind in the ABL. The 
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geostrophic forcing in turn affected the surface stress. The magnitude and variability 
of these interactions are important for understanding ice, ocean and atmospheric dy- 


namics.in the MIZ. 
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